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We have proposed the non-Abelian dual superconductivity picture for quark confinement in the 
SU(3) Yang-Mills (YM) theory, and have given numerical evidences for the restricted-field domi- 
nance and the non-Abelian magnetic monopole dominance in the string tension by applying a new 
formulation of the YM theory on a lattice. To establish the non-Abelian dual superconductivity 
| picture for quark confinement, we have observed the non-Abelian dual Meissner effect in the 5(7(3) 

, Yang-Mills theory by measuring the chromoelectric flux created by the quark-antiquark source, and 

the non-Abelian magnetic monopole currents induced around the flux. We conclude that the dual 
£S) ' superconductivity of the SU (3) Yang-Mills theory is strictly the type I and that this type of dual 

^ , superconductivity is reproduced by the restricted field and the non-Abelian magnetic monopole 

■ part, in sharp contrast to the SU(2) case: the border of type I and type II. 

Q \ PACS numbers: 12.38.Aw, 21.65.Qr 
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I. INTRODUCTION 

Oh 1 Quark confinement follows from the area law of the Wilson loop average. The dual superconductivity is the 
promising mechanism for quark confinement [l|. In many preceding works, the Abelian projection 0] was used to 

1 i ■ perform numerical analyses, which exhibited the remarkable results such as Abelian dominance [3] , magnetic monopole 

dominance [1, and center vortex dominance [6] in the string tension. However, these results are obtained only in 
T ~' . special gauges: the maximal Abelian (MA) gauge and the Laplacian Abelian gauge within the Abelian projection, 
Jft ' which breaks the gauge symmetry as well as color symmetry (global symmetry). 

In order to overcome the shortcomings of the Abelian projection, we have presented a new lattice formulation of 
ly— ^ ■ SUJN) Yang-Mills (YM) theory in the previous papers [7|, l|| (as a lattice version of the continuum formulations 
[9l.[l0| for SU(2) and [HI for SU(N)), which gives a decomposition of the gauge link variable suited for extracting the 
dominant modes for quark confinement in t he g auge independent way. In the case of SU(2), the decomposition of the 
gauge link variable was given on a lattice (12| - [l5j as a lattice version of the Cho-Duan-Ge-Faddeev-Niemi decomposition 

■ Q- For the gauge group G — SU(N) (N > 3), it was found that the extension of the decomposition from SU(2) to 
t-H ' SU(N) (N > 3) is not unique and that there are a number of possible ways of decompositions discriminated by the 

stability subgroup H of G, while there is the unique option of H = U(l) in the SU(2) case [16J. 

For the case of G = SU(3), in particular, there are two possibilities which we call the maximal option and the 
minimal option. The maximal option is obtained for the stability group H = U(l) x U(l), which enables us to give 
a gauge invariant version of the MA gauge as the Abelian projection |17l - [20| . The minimal one is obtained for the 
stability group H = U(2) = SU(2) x U(l), which is suited for representing the Wilson loop in the fundamental 
representation as derived from the non-Abelian Stokes theorem [2lL |22| . In the static potential for a pair of quark 
and antiquark in the fundamental representation, we have demonstrated in (23rl25| and [2fj| : (i) the restricted-field 
dominance or "Abelian" dominance (which is a gauge-independent (invariant) extension of the conventionally called 
Abelian dominance): the string tension cry obtained from the decomposed y-field (i.e., restricted field) reproduced the 
string tension (Tf u n of the original YM field, av/cfuii = 93 ± 16%, (ii) the gauge-independent non-Abelian magnetic 
monopole dominance: the string tension ay extracted from the restricted field was reproduced by only the (non- 
Abelian) magnetic monopole part tr mon , Cmon/ov = 94 ± 9%. 
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The purpose of this paper is to give further evidences for establishing the non-Abelian dual superconductivity 
picture for quark confinement in SU(3) Yang-Mills theory claimed in |26| by applying our new formulation to the 
SU(3) YM theory on a lattice. First, we study the dual Meissner effect by measuring the distribution of chromo-flux 
created by a pair of static quark and antiquark. We compare the chromo-flux of the original Yang-Mills field with that 
of the restricted field and examine if the restricted field corresponding to the stability group H = U(2) reproduces 
the dual Meissner effect, namely, the dominant part of the chromoelectric field strength of SU (3) Yang-Mills theory. 
Second, we measure the possible magnetic monopole current induced around the flux connecting a pair of static quark 
and antiquark. Third, we focus on the type of dual superconductivity, i.e., type I or type II. In the SU(2) case, 
the extracted field corresponding to the stability group H = U(l) reproduces the dual Meissner effect, which gives a 
gauge invariant version of MA gauge in the Abelian projection, as will be given in p7| . In this paper, we find that 
the dual superconductivity of the SU(3) Yang-Mills theory is indeed the type I, in sharp contrast to the SU(2) case: 
the border of type I and type II [27| . 

The chromo-flux of the Yang-Mills field has been already measured by using a gauge- invariant Wilson line/loop 
operator [28[ by the preceding works: [20M36I ] for SU(2) case and f35T - [56l | for the SU(S) case. However, most of the 
SU(2) case were done using the Abelian projection and there were no direct measurement of the dual Meissner effect 
in the gauge independent (invariant) way, except for quite recent studies [36|, HH, HI] • 



II. LATTICE FORMULATION 

We give a quick review of a new formulation of the lattice SU(N) YM theory 0, §]. In particular, the minimal 
option of the formulation enables us to extract the dominant mode of quark confinement, if we focus our attention 
on quarks belonging to the fundamental representation of the SU(N) gauge group |25|, Ha]- For the original SU(3) 
gauge link variable U x ^ G SU(3), we wish to decompose it into new variables V x M and X x M which have values in the 
517(3) group, i.e., e SU(3), V Xjfl G 5Z7(3): 

U X)li = X x ^V X4l G 5(7(3), (1) 

so that V x .^ could be the dominant mode for quark confinement, while X Xtfl is the remainder. In this decomposition, 
we require that V Xtfl is transformed in the same way as the original gauge link variable U Xifi and X x ^ as a site variable 
by the full SU(3) gauge transformation ft x : 

U x ,n > V x v = il x U Xifl il x+fl , (2a) 

v Xtfl — > V x , v = n^n^, x x ^ — > x' x>v = n x x x ^ x . (2b) 

First, we introduce the key variable called the color field. In the minimal option, the color field is defined by 

h x = C(A 8 /2)e f G Lie[SU(3)/U(2)], (3) 

with A 8 being the Gell-Mann matrix and £ the SU(3) group element. Once the color field is introduced, the above 
decomposition is obtained by solving the defining equation: 

Dl[V]h x := ~ [V^K+v - KV X J = 0. (4) 
This defining equation can be solved exactly Q, and the solution is given by 

x *,n = det^^) 1 /^" 1 , V x>il = Xl^U^ = g x L Xtll U Xflx , (5a) 

L x ,^ — (L X: ^L X ^) L x ,h, (5b) 



5 n 

l x ,v = -l + y -(K + u Xttl h. x+IJ ui tP ) + Nli,/; r .,,h,^ j( rj. 



(5c) 



Here the variable g x is the U(2) part undetermined from Eq.@ alone. In what follows, therefore, we put 

9x = 1, (6) 
so that the above defining equations Q and ^ correspond respectively to the continuum version (Tlj : 

iyv]h(a:) -0, (7a) 



X^x) - ~[h(x), [h(x),X^x)} = 0. (7b) 
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For the physical meaning of the defining equation, see [22]. In the naive continuum limit, indeed, ([5]) reproduces the 
decomposition in the continuum theory 



A^x) = V„(a;)+X M (x), 
4 
3 

X^z) = | [h(x), [h^A^fx)]] +i<T^ [d„h[.r). b(,-)] . 



V„U) = A„(.r) - - [h(.rl. hi.» i.A M (x)]]-iff- 1 | [d M h(x),h(x)] , 

,-i4 



(8b) 



Thus the decomposition is uniquely determined from Eqs.©, if the color field is specified. In order to determine 
the configuration {h x } of color fields, we use the reduction condition 0, H| which guarantees that the new theory 
written in terms of new variables {X Xttl ,V Xifl ,) is equipollent to the original YM theory. Here, we use the reduction 
condition: for a given configuration of the original link variables U Xl/J/ , color fields {h x } are obtained by minimizing 
the functional: 



F Ied [{h x }}=^{(D^[U x Jh x )HD^[U x Jh x )}. 



(9) 



From the non- Abelian Stokes theorem [2lj, [22| and the Hodge decomposition of the field strength [V] , we can 
define the magnetic part (magnetic monopole current) in the gauge independent way. The lattice version of the 
gauge-invariant magnetic-monopole current k is given by using the y-field as 



V x>IJ y x +wV£ +v li y£ v = exp (-ige 2 J' M1 ,[V](x)) , 



e 2 e^(x) := -arg Tr 



fe^(x) = 2im li (x) = d v *Q,j,u(x) := 



xpd v ® a p(x). 



(10a) 
(10b) 

(10c) 



Note that the the current k^ is the non- Abelian magnetic monopole current, since V x ^ is the decomposed restricted 
field. 



III. METHOD AND RESULTS 



We generate configurations of the YM gauge link variable {U x ^} using the standard Wilson action on a 24 4 lattice 
at (3 = 6.2. The gauge link decomposition is obtained according to the framework given in the previous section: the 
color field configuration {h x } is obtained by solving the reduction condition of minimizing the functional eq.© for 
each gauge configuration {U x _^\, and then the decomposed variables {V^ 1(tl }, {X XifJ } are obtained by using the formula 
eq.([5])- In the measurement of the Wilson loop average, we apply the APE smearing technique to reduce noises (57j . 



A. Dual Meissner effect 



We investigate the non- Abelian dual Meissner effect as the mechanism of quark confinement. In order to extract 
the chromo- field, we use a gauge-invariant correlation function proposed in |28| : The chromo-field created by a quark- 
antiquark pair in SU (N) Yang-Mills theory is measured by using a gauge-invariant connected correlator between a 
plaquette and the Wilson loop (see FigfT]): 

_ (tT(UptfWL)) 1 (tr(t/p)tr(WQ) 
PW (tr(WO) N (tr(VF)) ' { } 

where W is the Wilson loop in Z-T plane representing a pair of quark and antiquark, Up a plaquette variable as the 
probe operator to measure the chromo-field strength at the point P, and L the Wilson line connecting the source W 
and the probe Up. Here L is necessary to guarantee the gauge invariance of the correlator pw and hence the probe is 
identified with LUpL^. The symbol (O) denotes the average of the operator O in the space and the ensemble of the 
configurations. In the naive continuum limit e — ► 0, indeed, pw reduces to the field strength in the presence of the qq 
source: 

2 (tr(ige*L^W)) 



(tr(W0) 
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FIG. 1: (Left) The setup of measuring the chromo-fiux produced by a quark-antiquark pair. (Right) The gauge-invariant 
connected correlator (UpLWL^) between a plaquette U and the Wilson loop W . 
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color flux: restricted field (L/e =8 , z/e =4 ) 
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FIG. 2: Measurement of components of the chromoelectric field E and chromomagnetic field B as functions of the distance y 
from the z axis. (Left panel) the original SU(3) YM field, (Right panel) the restricted U(2) field. 



where we have used U Xtfl = exp(— igtA^(x)) and hence Up = exp(— igt 2 J-^ v ). Thus, the gauge-invariant chromo-field 
strength F^ v [U] is given by 



F^[U] :=e 




Pw, 



(13) 



where j3 := 2N/g 2 is the lattice gauge coupling constant. Note that the connected correlator pw is sensitive to the 
field strength, while the disconnected one probes the squared field strength: 



Pw : = 



(tr (W) tr (U p )) 
(tr(W)) 



(tr(U p )) 



/JF2 \ _/j2\ 



(14) 



We measure correlators between the plaquette Up and the chromo-field strength of the restricted field V x ^ as well 
as the original YM field U X4l . See the left panel of Fig. [1] Here the quark and antiquark source is introduced as 8 x 8 
Wilson loop (W) in the Z-T plane, and the probe (U p ) is set at the center of the Wilson loop and moved along the 
Indirection. The left and right panel of Fig. [5] show respectively the results of measurements for the chromoelectric 
and chromomagnetic fields F^lU] for the original SU(3) field U and F^IV] for the restricted field V, where the field 
strength F MI ,[V] is obtained by using V x ^ in eq(JTTj) instead of U x ^\ 



(15) 
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FIG. 3: The distribution in Y-Z plane of the chromoelectric field E z connecting a pair of quark and antiquark: (Left panel) 
chromoelectric field produced from the original YM field, (Right panel) chromoelectric field produced from the restricted U(2) 



We have checked that even if W[U] is replaced by W[V], together with replacement of the probe LUpL^ by the 
corresponding V version, the change in the magnitude of the field strength remains within at most a few %. 

From FigJ5]we find that only the E z component of the chromoelectric field (E x , E y , E z ) — (Fiq, F20, F30) connecting 
q and q has non-zero value for both the restricted field V and the original YM field U. The other components are zero 
consistently within the numerical errors. This means that the chromomagnetic field (B x , B y , B z ) — (F23, -F31, -F12) 
connecting q and q does not exist and that the chromoelectric field is parallel to the z axis on which quark and 
antiquark are located. The magnitude E z quickly decreases in the distance y away from the Wilson loop, as discussed 
in more detail in subsection C. 

To see the profile of the nonvanishing component E z of the chromoelectric field in detail, we explore the distribution 
of chromoelectric field on the 2-dimensional plane. Fig.[3]shows the distribution of E z component of the chromoelectric 
field, where the quark-antiquark source represented as 9 x 11 Wilson loop W is placed at (Y, Z) = (0, 0), (0, 9), and the 
probe U is displaced on the Y-Z plane at the midpoint of the T-direction. The position of a quark and an antiquark 
is marked by the solid (blue) box. The magnitude of E z is shown by the height of the 3D plot and also the contour 
plot in the bottom plane. The left panel of Fig. [3] shows the plot of E z for the SU (3) YM field U, and the right panel 
of Fig. [3] for the restricted-field V. We find that the magnitude E z is quite uniform for the restricted part V, while it 
is almost uniform for the original part U except for the neighborhoods of the locations of q, q source. This difference 
is due to the contributions from the remaining part X which affects only the short distance, as will be discussed in 
the next section. 



Next, we investigate the relation between the chromoelectric flux and the magnetic current. The magnetic(- 
monopole) current can be calculated as 



where F[V} is the field strength (2-form) of the restricted field (1-form) V, d the exterior derivative and * denotes the 
Hodge dual operation. Note that non-zero magnetic current follows from violation of the Bianchi identity (If the field 
strength was given by the exterior derivative of V field (one-form), F\V] = dV, we would obtain k =*c? 2 V = 0). 

Fig. [4] shows the magnetic current measured in X-Y plane at the midpoint of quark and antiquark pair in the 
Z-direction. The left panel of Fig. 2] shows the positional relationship between chromoelectric flux and magnetic 
current. The right panel of Fig. @] shows the magnitude of the chromoelectric field E z (left scale) and the magnetic 
current k (right scale) . The existence of nonvanishing magnetic current k around the chromoelectric field E z supports 
the dual picture of the ordinary superconductor exhibiting the electric current J around the magnetic field B. 

In our formulation, it is possible to define a gauge-invariant magnetic-monopole current fc M by using I^-field, which 
is obtained from the field strength J-^fV] of the field V, as suggested from the non-Abelian Stokes theorem [2ll. I22I]. 
It should be also noticed that this magnetic-monopole current is a non-Abelian magnetic monopole extracted from 
the V field, which corresponds to the stability group H = U(2). The magnetic-monopole current fc M defined in this 



field. 



B. Magnetic current 



k = *dF[V] 



(16) 



6 




0.07 

0.06 

0.05 

0.04 
Eze 2 
0.03 

0.02 

0.01 





' ' chrom'o flux Ez' +~ 

magnetic current Jm x 



x + 



+ x 



t >< + X ± v 



4 6 8 

distance y/e 



10 



12 



0.035 
0.03 
0.025 
0.02 

Jme 3 
0.015 

0.01 

0.005 





FIG. 4: The magnetic-monopole current k induced around the flux along the z axis connecting a quark-antiquark pair. (Left 
panel) The positional relationship between the chromoelectric field E z and the magnetic current k. (Right panel) The magnitude 
of the chromo-electronic current E z and the magnetic current J m — |k| as functions of the distance y from the z axis. 

way can be used to study the magnetic current around the chromoelectric flux tube, instead of the above definition 
of k p6p . The comparison of two monopole currents k will be done in the forthcoming paper. 



C. Type of dual superconductivity 

Moreover, we investigate the QCD vacuum, i.e., type of the dual superconductor. The left panel of Figf5]is the plot 
for the chromoelectric field E z as a function of the distance y in units of the lattice spacing e for the original SU (3) 
field and for the restricted field. 

In order to examine the type of the dual superconductivity, we apply the formula for the magnetic field derived 
by Clem (59} in the ordinary superconductor based on the Ginzburg-Landau (GL) theory to the chromoelectric field 
in the dual superconductor. In the GL theory, the gauge field A and the scalar field cf> obey simultaneously the GL 
equation: 



(17) 



and the Ampere equation: 



d'F^ + *g[0*(0 M - iqA^) - (8^ - iqA^)*4>] = 0. (18) 

Usually, in the dual superconductor of the type II, it is justified to use the asymptotic form Ko(y/X) to fit the 
chromoelectric field in the large y region (as the solution of the Ampere equation in the dual GL theory). However, it 
is clear that this solution cannot be applied to the small y region, as is easily seen from the fact that K (y/X) — > oo 
as y — > 0. In order to see the difference between type I and type II, it is crucial to see the relatively small y region. 
Therefore, such a simple form cannot be used to detect the type I dual superconductor. However, this important 
aspect was ignored in the preceding studies except for a work [361 ] . 

On the other hand, Clem [5{| does not obtain the analytical solution of the GL equation explicitly and use an 
approximated form for the scalar field <fi (given below in (|20jl ). This form is used to solve the Ampere equation 
exactly to obtain the analytical form for the gauge field A^ and the resulting magnetic field B. This method does not 
change the behavior of the gauge field in the long distance, but it gives a finite value for the gauge field even at the 
origin. Therefore, we can obtain the formula which is valid for any distance (core radius) y from the axis connecting 
q and q: the profile of chromoelectric field in the dual superconductor is obtained: 



(19) 
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FIG. 5: (Left panel) The plot of the chromoelectric field E z versus the distance y in units of the lattice spacing e and the 
fitting as a function E z {y) of y according to (|23[) . The red cross for the original SU(3) field and the green square symbol for 
the restricted field. (Right panel) The order parameter (f> reproduced as a function cf>(y) of y according to (|2Up . togather with 
the chromoelectric field E z (y). 
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1.828 ±0.023 


3.26 ±0.13 


3.84 ±0.19 


2.96 ±0.3 
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TABLE I: The properties of the Yang-Mills vacuum as the dual superconductor obtained by fitting the data of chromoelectric 
field with the prediction of the dual Ginzburg-Landau theory. 



provided that the scalar field is given by (See the right panel of FigJS]) 

<P(y) = — -4= ; V (20) 

where K v is the modified Bessel function of the z/-th order, A the parameter corresponding to the London penetration 
length, £ a variational parameter for the core radius, and $ external electric flux. In the dual superconductor, we 
define the GL parameter k as the ratio of the London penetration length A and the coherence length £ which measures 
the coherence of the magnetic monopole condensate (the dual version of the Cooper pair condensate): 

(21) 

It is given by [H^] 

/ t = V2^ v /l-i < : 2 (C/A)/^ 1 2 (C/A). (22) 

According to the formula Eq. (JTOJ) , we estimate the GL parameter n for the dual superconductor of SU(3) YM 
theory, although this formula is obtained for the ordinary superconductor of U(l) gauge field. By using the fitting 
function: 

$ 1 1 ,1 C 



we obtain the result shown in Table HI The superconductor is type I if K < K c , while type II if K > K c , where the 
critical value of GL parameter dividing the type of the superconductor is given by k c — l/v2 ~ 0.707. Our data 
clearly shows that the dual superconductor of SU(3) Yang-Mills theory is type I with 

k = 0.45 ±0.01. (24) 

This result is consistent with a quite recent result obtained independently by Cea, Cosmai and Papa [36j]. Moreover, 
our result shows that the restricted part plays the dominant role in determining the type of the non-Abelian dual 
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FIG. 6: The rescaled correlation correlation functions r 3/2 (O(r)O(O)) for O = A, V, X for 24 4 lattice with /3 = 5.7, 5.85, 6.0. 
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The physical scale is set in units of the string tension 0" p h ys . The correlation functions have the profile of cosh type because of 
the periodic boundary condition, and hence we use data within distance of the half size of lattice. 

superconductivity of the SU(3) Yang-Mills theory, reproducing the same result jH^JHH, i.e., type I with 

k = 0.48 ±0.02. (25) 

This is a novel feature overlooked in the preceding studies. Thus the restricted-field dominance can be seen also 
in the determination of the type of dual superconductivity where the discrepancy is just the normalization of the 
chromoelectric field at the core y — 0, coming from the difference of the total flux $. These are gauge-invariant 
results. Note again that this restricted-field and the non-Abelian magnetic monopole extracted from it reproduce the 
string tension in the static quark-antiquark potential [25|, [2(| . 

Our result should be compared with the result obtained by using the Abelian projection: Y.Matsubara et. al [35j ] 
suggests k = 0.5 ~ I(which is j3 dependent), border of type I and type II for both SU{2) and SU(3). In SU(2) case, 
on the other hand, there are other works [33, H3| which conclude that the type of vacuum is at the border of type I 
and type II. 

We should mention the work [43| which concludes that the dual superconductivity of SU(3) Yang-Mills theory is 
type II with k = 1.2 ~ 1.3. This conclusion seems to contradict our result for SU(3). If the above formula (fT9|) is 
applied to the data of [43J, we have the same conclusion, namely, the type I with k = 0.47 ~ 0.50. Therefore, the 
data obtained in [43[ are consistent with ours. The difference between type I and type II is attributed to the way of 
fitting the data with the formula for the chromo-field. 



IV. GLUON PROPAGATORS 



Finally, we study the 2-point correlation functions (propagators) of the new variables and the original Yang-Mills 
field variables, which are defined by 

D 00 (x~y) := (OA(x)0*(y)) for 0£(x) € {V£ v , X^, A*J, (26) 

where an operator O^x) is defined by the linear type, e.g., A x >^ :— {U x ^ — U^. ^)/(2ige) where x' means the mid- 
point of x and x + efi. In order to calculate the propagators, we must impose a gauge fixing condition, and we have 
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adopted the lattice Landau gauge (LLG). Fig. [6] shows the logarithmic plot of the scaled propagators r 3 / 2 Doo(r) 
as a function of r = \x — y\, where the distance r is drawn in units of the string tension (Xphys, and data of lattice 
spacing is taken from the TABLE I in Ref.[58]. The propagator Dyy damps slowly and has almost the same damping 
behavior as Daa, while the Dxx damps quickly. Thus, from the viewpoint of the propagator, the y-field plays the 
dominant role in the deep infrared region or the long distance, while X-field is negligible in the long distance. The 
rapid disappearance of X contribution in the long distance is helpful to understand the difference of the profile of the 
flux tube in Fig|3] 

We estimate the mass Mo (i.e., the rate of exponential fall-off) from the propagator Doo( r ) by using the Fourier 
transformation of the massive propagator in the Euclidean space, which behaves for large Mor as 

^-f&^j^*^!^ (Mor>>1) - (27) 

and hence the scaled propagator r 3 ^ 2 Doo ( r ) should be proportional to exp(— Mor) . In order to perform the parameter 
fitting of Mo for O = {~V x ',fj,, A x / jAl }, we use data in the region [2.0, 4.5] and exclude the data near the midpoint of the 
lattice to eliminate the finite volume effect, while for O — X x ' jjU we use the region [1.0, 3.5]. We use CTphys = MOMeV 
to obtain preliminary result: 

M A ~ 0.76 GeV, M v ~ 0.73 GeV , M x ^ 1.15 GeV, (28) 

which should be compared with result of the maximal option [171 ] in LLG, and also result of the Abelian projection 
in the maximal Abelian gauge p~8( ] . 

V. SUMMARY AND OUTLOOK 

In this paper, we have given further numerical evidences for confirming the non-Abelian dual superconductivity 
for SU(3) YM theory proposed in (2(|. For this purpose, we have used our new formulation of SU(3) YM theory 
on a lattice 0, H| to extract the restricted field from the original SU(3) YM field, which has played a dominant role 
in confinement of quarks in the fundamental representation, i.e., the restricted-field dominance and the non-Abelian 
magnetic monopole dominance in the string tension, as shown in the previous studies |26j . 

We have focused on the dual Meissncr effect and have measured the chromoelectric field connecting a quark and an 
antiquark for both the original YM field and the restricted field. We have observed the dual Meissner effect in SU(3) 
YM theory, i.e., only the chromoelectric field exists and the magnetic- monopole current is induced around the flux 
connecting a quark and an antiquark. Moreover, we have determined the type of non-Abelian dual superconductivity, 
i.e., type I for SU(3) YM theory, which should be compared with the border of type I and II for the SU(2) YM 
theory. These features are reproduced only from the restricted part. These results confirm the non-Abelian dual 
superconductivity picture for quark confinement. 

In order to confirm the type of the non-Abelian dual superconductivity, it will be interesting to investigate if the 
interactions between two chromoelectric flux tubes are attractive (expected for the type I) or repulsive (expected for 
the type II). The interactions should reflect the internal non-Abelian nature which depends on the distance between 
two tubes. It will be also interesting to extend our result to the finite temperature case to see the fate of the dual 
superconductivity, and to see when and how the dual superconductivity disappears above the critical temperature. 
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